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Magnesium diboride is remarkable not only for its unusually high transition 
temperature of 39K for a sp-bonded metal [|l| and its reduced isotope-effect ex- 
ponent but also for its unusual superconducting gap behavior [ ^16|| . Some 
earlier experiments seem consistent with conventional single-gap superconduc- 
tivity but recent measurements indicate the existence of more than one 
gap |9|- |16|l - For example, specific heat measurements show an anomalous large 
bump at low temperature, inconsistent with single-gap BCS theory. Moreover, 
there is a significant variation in the measured values for the gaps. Here, we re- 
port first-principles calculations of the k- and T-dependent superconducting gap 
A{k,T) in MgB2 and its manifestation in various measured quantities. Because 
its Fermi surface has disconnected sheets with very diff"erent electron-phonon 
coupling strengths, our calculations show that near T = 0, the values of A{k) 
cluster into two groups: large values (~ 6.5 to 7.5 meV) on the strongly coupled 
sheets and small values (~ 1 to 3 meV) on the weakly coupled sheets. The 
calculated gap, quasiparticle density of states, and specific heat and their tem- 
perature dependences are in agreement with the recent measurements which 
support that MgB2 is a multiple gap superconductor. In fact, theory predicts 
four prominent values for the gap at low T. 

Most theoretical work on MgB2 has been focused on the electronic structure, the phonon 



structure, and the electron-phonon interaction |T7|-p3|. It is shown that the electron-phonon 
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interaction in MgB2 varies strongly on the Fermi surface pT|- p^] . Also, recent two-band 



model 1^ and first-principles ||2^ studies of the transition temperature support a multi- 



gap scenario However, no detailed quantitative calculations have been presented for 

the temperature and fc-space dependence of the superconducting gap. The calculational 
framework P3 used here is based on the strong coupling formalism of superconductivity 



established by Eliashberg |^5|,^. The Eliashberg formalism is a more general case of the 



original formulation of the BCS theory for phonon-mediated pairing, taking into account 
strong electron-phonon coupling. When applying the strong coupling theory to MgB2, we 
include the momentum dependency of the anisotropic Eliashberg equation and also take 
into account the anharmonic effect on phonon frequencies. Material-dependent parameters 
in the equation are obtained by ab initio pseudopotential density functional calculations. 
This approach has already been shown to describe the superconducting transition of MgB2 
very accurately; the obtained specific heat mass enhancement factor A = 0.61, the transition 
temperature = 39K, and isotope-effect exponent ob = 0.31 are all in excellent agreement 



with experiments Here, we compute the gap function A(k,uj) and the renormalization 
factor Z{k,uj) from the fully anisotropic Eliashberg equation at various temperatures. This 
allows us to address all the superconducting properties except those related to the presence 
of an applied magnetic field. 

Figure 1 shows the calculated superconducting gap A(fc) on the Fermi surface at 4 K. 
The Fermi surface of MgB2 consists of four sheets: two-dimensional (2D) light hole and 
heavy hole sheets forming coaxial cylinders along T to A, a three-dimensional (3D) hole 
sheet connecting regions near K and M, and a 3D electron sheet connecting regions near 
H and L. The 2D light and heavy hole sheets are derived from cr-antibonding states of 
boron p^^y orbitals, while the 3D hole and electron sheets are derived from vr-bonding and 
antibonding states of boron orbitals, respectively. The calculated density of states at the 
Fermi energy is 0.115 states/eV- atom- spin; 44 % of which comes from the 2D a cylindrical 
sheets and the rest comes from the 3D tt sheets. The superconducting gap is nonzero 
everywhere on the Fermi surface. Thus, the symmetry of the superconducting gap is s- 



wave, but the size of the gap changes on the Fermi surface. The gap values cluster into two 
groups. The largest gap is on the 2D light hole a cylindrical sheet (shown red in Fig. 1), 
where the average gap is 7.2 meV with variations of less than 0.1 meV. On the 2D heavy hole 
a cylindrical sheet (shown orange in Fig. 1), the superconducting gap ranges from 6.4 meV 
to 6.8 meV, having an average of 6.6 meV, with maximum value near F and minimum value 
near A. The average of the gap values on the two 2D a cylindrical sheets is 6.8 meV. The 
superconducting gap is significantly smaller and more spread out on the 3D vr sheets (shown 
green and blue in Fig. 1), ranging from 1.2 meV to 3.7 meV. Out of the two 3D vr sheets, 
the hole sheet shows a slightly larger gap with a peak at 2.1 meV and the electron sheet 
shows a somewhat smaller gap with a peak at 1.5 meV. The average of the gap values on the 
two 3D vr sheets is 1.8 meV. This very large variation in the gap value originates from the 
differences in electron-phonon coupling strengths on the different sheets of the Fermi surface. 
Our ah initio calculations p3| demonstrate that the electron-phonon coupling strength is 4 



to 5 times stronger on the cylindrical a sheets (exceeding 2.5 for some states) than on the 
vr sheets. Our result is consistent with the recent experiments reporting two gaps [^-[T^ 
because the superconducting gap in MgB2 can be grossly viewed as consisting of a large gap 
of ~ 6.8 meV on the strongly coupled 2D cr-cylindrical sheets and a small gap of ~ 1.8 meV 
on the weakly coupled 3D vr sheets. However, this simplication is limited. As seen in Fig. 1, 
the large and small gaps both have substantial variation in value: 6.4 to 7.2 meV on the 2D 
(J sheets and 1.2 to 3.7 meV on the 3D vr sheets. The experimental interpreted gaps have 
a range of 1.5 to 3.5 meV for the small gap and a range of 5.5 to 8 meV for the large gap 



Figure 2 shows the temperature dependence of the superconducting energy gap on the 
Fermi surface from 4K to 38K which may be probed in tunneling, optical, and specific heat 
measurements. The vertical blue curves in Fig. 2 present the distribution of the super- 
conducting gap values at various temperatures, and the red lines are curves of the form 



A(T) = A(0)y 1 — {T/Tc)P fitted separately to the averaged values of the calculated gap on 
the 2D a sheets and those on the 3D vr sheets. Our calculated results give A(0) = 6.8 meV 



{2A{0) /keTc = 4.0) for the averaged gap at T = on the 2D a sheets, and A(0) = 1.8 meV 
{2A{0) /ksTc = 1.06) for that on the 3D vr sheets. The theory shows that the superconduct- 
ing gap opens at Tc throughout the Fermi surface, and that the gap widens much faster on 
the 2D a sheets than on the 3D vr sheets at T ^ Tc. 

From the calculated gap function A{k,uj), the quasiparticle density of states is given by 

N{u;)/N{0) = Re ( , ^ , ) , (1) 

\^J{uj + iT)^ -A{k,uj)y 

where (■ ■ ■) indicates an average over a surface of constant u). Figure 3 depicts the theoret- 
ical quasiparticle density of states calculated with an assumed finite lifetime F of 0.1 meV. 
The quasiparticle density of states at 4K shows three discernable peaks. One peak is at 2.2 
meV, which is slightly greater than the dominant size of the gap (2.1 meV) on the 3D vr 
hole sheet. The other two peaks are at 6.8 and 7.2 meV, which are the maximal sizes of the 
superconducting gap on the 2D heavy and light hole vr sheets, respectively. The quasiparticle 
density of states can be deduced from tunneling experiments and various spectroscopic mea- 
surements [p HTTl , p!5| , p!6| , but a direct quantitative comparison requires knowledge of various 
physical parameters involved in a specific experiment. 

A fundamental measurement which can provide a probe of the superconducting quasi- 
particle spectrum is the temperature-dependent specific heat. Figure 4 shows the calcu- 
lated T-dependent specific heat of MgB2 through the superconducting transition. For the 



normal-state specific heat Cn = 7T, we obtain 7 = 2.62 mJ/mol-K p3[, which agrees 



well with experimental values 2.6 mJ/mol-K^ |jT3| and 2.7 mJ/mol-K^ [1T2| , [T^ . For the 



superconducting-state specific heat Cs, we first calculate the free energy difference (F5 — Fj 



N) 



of the superconducting and normal states pSf and then obtain the specific heat difference 
by 

Cs-Cm = -T^^{Fs-F^) . (2) 



The measured low T specific heat ||T^^H[ shows substantial magnitude and a large bump at 



about lOK. This anomalous behavior is inconsistent with a 1-gap BCS model. However, our 



results are in excellent agreement with experiment. The origin of the bump in our calculated 
curve is the existence of low energy excitations above the small gap on the weakly coupled 
3D TT sheets of the Fermi surface. The overall shape and magnitude of the calculated specific 
heat curve agrees very well with experiments, especially below 30K, and shows a jump of 
AC/7TC = 1.0 at Tc which is within the range of the measured values. 
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FIG. 1. (color) The superconducting energy gap on the Fermi surface at 4 K. Data is given in 
a color scale (top panel). The bottom panel depicts the distribution of gap values and the color 
scale used. 
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FIG. 3. The quasiparticle density of states at various temperatures. 
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